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Abstract: The fast direct solver based on skeletonization is extended to the method of moments 
(MoM) solution of volume integral equations (VIE). The impedance matrix is approximated by the 
hierarchical data sparse representation through skeletonization. The inverse matrix is then implicitly 
computed. Numerical experiments are conducted to study the performance of the direct solver. 
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I. Formulation of VIE 


ouppose an inhomogeneous, linear and — 
isotropic object with volume V in free EC 


space, the total electric field E with the 
time factor of e" can be expressed 
asl!!: 


E = E" +L Ji 1o = y Lo / Eo 


L{x}(r)= - jk, | av L1 es VV ponatur) 





I 
e(r) 


Discretization and testing 


Meshing the object by tetrahedrons, the unknown flux D can be represented by 
the SWG basis function [! . And with the Galerkin scheme, the resultant linear 


system reads: 
D(r)- »,, D,f,(r) 


N The number of unknowns 








f „ The n-th SWG basis function 


D Associated unknowns coefficient 
£ The resultant impedance matrix Z is fully populated 
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II. RS-Based Direct Solver 


[he skeletonization framework 


Suppose the field group o and source group s, contains N, testing RWG 
functions and M, basis RWG functions. The associated coupling matrix between 
these two groups Z,, is a low rank matrix . 


Z with dimension N, x N Skeletonization 


sain fog Sho So R 
Decomposition (ID)[1] OS O OS OS ID for left decomposition 
° finding principle 


C — E kxk, <N, 
basis/testing functions; /. M 


conducting under the (ds. S TOW PIRRE matrix 

control of error ‘of y" , with receiving 

threshold ME skeleton HWGs in kK, 
* recursive decomposition group o. 


would make use of rank L pm matrix from S 
deficiency to the most. to S. 














Š with dimension N, x K, 





X "NN Sksletonized 
" matrix of Z,, , with K, 
radiation skeleton RWGs 
In group S. 
R. projection matrix from S ` 
to Z... 
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[he skeletonization framework 








1,8 
7 7 mxn| [Low rank block 
Mj D 
[== i 
: : i Skeletonization via ID 
Zi D, Uga Ti 66 2 AR SU k.«n. 
[] l [] J i i J J 


Z 
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Inverse of the impedance matrix 


The inverse of the impedance matrix 





N number of RWG basis functions 
K total rank of all groups 


D,E, F,G are all block diagonal matrices 


To systematically take advantage of the ran deficiency of the couplings, a 
hierarchical decompostition of the original geometry would be desirable. 


Z'=E(S+D) F+G > z^ -E-(E, (s,+B,) -F, t6) F +G 





lo sum up, the | RIT: Matrix Direct 


Skeleton Inversion Solver 


whole process Is 
In two stages: 





III. The strategy of selecting 
skeletons 


To figure out the skeleton in 
group i for all receiving groups, a 
proxy matrix Z’ should be 
constructed!!! straightforwardly 
from Z. Despite high accuracy, It is 
of little practical use due to the 
O(N*) complexity of filling Z. A 
remedy is the Huygens principle 
|1]. [2], It replaces the observation Bisdhi 
points exterior to the objective Exterior Exterior 
group with a set of proxy points on 
an artificial cubic or spherical shell 
enclosing this group. (b) 
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Spherical shell 





Cubic shell 


Group Group 








IV. Numerical experiments 


Mesh Density 





N 11 010 2529) 44309 62 900 129211 
Example |: fixed frequency and D” (1.0 m), decreasing h 

I. 7 270 9 6/0 12 58 14 724 22 980 26 143 

Ts 34 | 05 254 538 2 /49 6 122 

Ti 195 415 954 I 574 7 215 10 476 

Mem. 622 | 220 2 263 3 26] 8 /34 12 687 

Example 2: fixed frequency and h (& 1.2 x 1077Agq), increasin E 


o 


Š XL 


200 639 

















n 5 634 8 636 12 040 13 804 20 895 21 020 


Ts 34 | 05 250 519 2 851 6 217 
Tit | 85 39] 903 | 462 7 OOS 9 023 
Mem. 450 1 021 2 126 3 149 7 907 9 792 


DO: target size Mc(I.e. h): mesh size 





Example 1 
fixed DO: 





Example 2 


fixed Mc: 








Time (s) or Memory (MB) 





—#— Time, fixed DO —t—- Mem., fixed DO 

—— Time, fixed Mc —o— Mem., fixed Mc 

Ideal O(N) — — - |deal O(N^1.3) 
*—*|deal O(N^1.5) 


Bi-static RCS and electric field: 


4 | Electric field in the X-Y plane 


8... Permittivity e 


0.90 














19 2.0 
m 0.84 Diameter 
— ° Z 104 
£ C—Á— C — m... F pin 
9 MIE Fast solver 7. 
© ó, <1.2x10 2 2 
"x 5. 
m 6, <1.0x10™ 
-10 





Permittivity 


—A— Time, total —o-— Mem. 
Ideal O(N) Ideal O(N^1.3) 
Fixed frequency (D°=0.1 A ) 
Fixed mesh 1.2x10 4, 
Decreasing £, 
Permittivity €, 


5.0, 8.0, 12.0, 15.0, „á 
21.0,25.0,30.0 -7% 





Error 


ó, «899x107 
0. <4.2x10~ 


Time(s) or Memory(MB) 


Bi-static RCS and electric field: 


Bi-static RCS(dB) 


-30 


-10 hace: Citi, 









-50 
-60 
-/0 : 
Fast solver ? e —30.0, MIE 
T - i e =30.0, Fast Solver 
EE : NN t |--..... e =5.0, MIE 
0.076 0.17 0.22 





Electric field in the X-Y plañe | 579-9. Fast Solver 
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Metallic nano-particles ^ —(—66.219,5.7015) 
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Bi-static RCS(dB) 


Fast solver 








aE ë 
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^9 -| Electric field in the X-Y plane * 
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Metallic nano-particles E, = (—66.219,5.7015) 





-47.0 
ann Tst: 1296s 
47.24 4 - % | Mem.: 3674MB 
cn x 
F -47.3 2 63/86 SWGs | 
O Š 15368 
F 47.4 skeletons 
m 
o 475 
Q 
C 
© 
= -476 
-47.7 
B 
0.023 0,038 0,059 0.071 
-47.8 
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V. Conclusions 


» [he direct solver based on skeletonization 
can be generally divided into two stages: 
selecting skeletons and computing the 
inverse matrix. 





» [he performance of the fast direct solver is 
investigated in terms of near field as well 
as far field computation. 


> Numerical experiments indicate that the 
proposed fast direct solver is efficient and 
accurate. 
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